FULL RANGE OF BLOW UP EXPONENTS FOR THE QUINTIC 
WAVE EQUATION IN THREE DIMENSIONS 



JOACHIM KRIEGER, WILHELM SCHLAG 

Abstract. For the critical focusing wave equation nu = u 5 on R 3+1 in 
the radial case, we prove the existence of type II blow up solutions with 
scaling parameter A(f) = f -1-1 ' for all v > 0. This extends the previous 
work by the authors and Tataru where the condition v > j had been 
imposed, and gives the optimal range of polynomial blow up rates in 
light of recent work by Duyckaerts, Kenig and Merle. 



1. Introduction 

We consider the energy critical focussing wave equation 

nu = u 5 ,n = d] - A (1.1) 

on R 3+1 , in the radial case. This equation has been intensely studied in a 
number of recent works: the remarkable series of papers (3) - J6J established 
a complete classification of all possible type II blow up dynamics, without 
proving their existence. In the works JU, 12) a constructive approach to 
actually exhibit and thereby prove the existence of such type II dynamics 
was undertaken. Recall that a type II blow up solution u(t, x) with blow up 
time T* is one for which 

limsup (\\u(t, + \\Ut(t, -)Il 2 ) < 00 

In |6], it is demonstrated that such solutions can be described as a sum of 
dynamically re-scaled ground states 



+W(x) = +(1 



2, _1 

2 



3 

plus a radiation term. In particular, for solutions where only one such bulk 
term is present, one can write the solution as 

u(t,x) = W A{t) (x) + w(t,x) + o Hl (l), W A = A5W(Ax), w{t,x) e H 1 (1.2) 
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where the "error" (w(t, •), dtw(t, •)) converges strongly in H 1 x L 2 as t — > T*, 
and we have the dynamic condition 

lim (T* - t)A{t) = oo (1.3) 

i — > 7 7 * 

In |8|, it was shown that such solutions with A(f) = f~ 1_v do exist, where 
v > A is arbitrary. This left the question whether for polynomial rates the 
condition < |1.3t is indeed optimal. Here we show that it is. 

Theorem 1.1. Let v > be given. Then there exists an energy class solution 
u(t,x), which in fact has regularity H 1+ i~, of the form (|1.2|), with 

A(£) = r 1_v 

Our method of proof is closely modeled on the construction from (8J, of 
which we now recall the main steps: 

(i) We write Uo(t, r) = W^( t ) (r) and iteratively modify Uo in the form 

2k-l 

«2fc-l = "o + Yj v i 
;=i 

so that «2)c-i satisfies (|l.lj) up to an error of size t N as f — > 0+; here 
N can be made as large as desired by taking A; large, and the size 
is measured relative to the energy inside a light cone with tip at 
r = 0,£ = 0. 

(ii) We seek an exact solution via a perturbation: u = «2ic-i + £- To solve 
for £ we switch to coordinates R = A(t)r, z = $ t A(s) ds = \t~ v . The 
variable t varies in the range To ^ t < oo. 

(iii) In the new coordinates, the driving linear operator is 

£ = -d RR - 5W 4 (R) on (0,00) 

We perform a spectral analysis of the operator, which exhibits a 
unique and simple negative eigenvalue, as well as continuum spec- 
trum; in addition, there is a zero energy resonance. The latter 
renders the spectral measure singular at zero energy. 

(iv) A contraction argument is set up for e with a vanishing condition 
at t = oo. For the contraction it is important that N in the first step 
is sufficiently large. 

The first three steps in this paper are essentially the same as in (8J. It 
is in the final step that we improve on the procedure in |5I. In fact, in 
Proposition 2.8 of that paper the condition v > \ arises in order to make the 
embedding 

(H 1+2 *(R 3 )) 5 c H 2a (R 3 ) 

for v/2 > la ^ |. This was used to control the quintic terms in the 
construction of the exact solution via iteration and application of a suitable 
parametrix. In fact, the difference shall consist in a more detailed analysis of 
the first iterate for the exact solution, which we exhibit as a sum of two terms, 
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one of which is smooth, the other of which satisfies a good L°°-bound near the 
origin. This latter feature comes from the fact that the loss of smoothness of 
the approximate solution occurs precisely on the characteristic light cone, 
and thus one expects the exact solution to be smoother near the spatial 
origin. 

2. The approximate solution for a power-law rescaling 

2.1. Generalities. In radial coordinates, <(!.!)> becomes 

2 

■LquinticU := U tt ~ U„ - -U r - U 5 = (2.1) 

This equation is known to be locally well-posed in the space 'H : = H 1 x 
L 2 (R 3 ), meaning that if (w(0), u t (0)) e ( H, then there exists a solution locally 
in time and continuous in time taking values in *H. Solutions need to be 
interpreted in the Duhamel sense: 

«m - cos (t |vi,/ + + [ ^«'-^i) „ 5(s) ds M 

These solutions -Cq U intk{ u ) = have finite energy: 

f 1 i/ 6 
E(u,Ut) = [^( w ^ + |Vu| 2 ) — — ] dx = const 
Jrs 2 6 

A special stationary solution is W(r) = (1 + r 2 /3)~5. By scaling, ?\h W(Ar) is 
also a solution for any A > 0. We are interested in letting A depend on time. 
More precisely, we would like to find solutions -Cq U intk u = of the form 

u{t,r) = W A(t) (r) + w{t,r), A(f) — oo as t 0+ (2.3) 

and w small in a suitable sense. It suffices to show that w remains small in 
energy, since this ensures that the solution blows up at time t = by the 
mechanism of "energy concentration" at the tip of the light-cone (f, r) = 
(0, 0) (think of solving backwards in time). 

2.2. The bulk term. For the convenience of the reader, and in order to 
correct some minor inaccuracies in |[8) such as the omission of harmless 
logR factors, we redo the iterative construction from that paper which 
constitutes step (i) from the four step procedure outlined above. We will 
henceforth fix A(f) = t~ l ~ v with v > 0. Set 

UQ (t,r) = A(t)*W(rA(*)) = A(t)*W(R) (2.4) 

While i/o is very far from being an approximate solution, the authors to- 
gether with D. Tataru showed in (81 that one can add successive corrections 

u k = UQ + v 1 + v 2 +v 3 + ... + v k (2.5) 

so that this function approximately solves (|2.1|) . To be specific, they achieved 
that ~CquintkU(t) goes to zero like f N in the energy norm restricted to a light 
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cone where N can be made arbitrarily large by taking k large. This is an 
iterative construction. Moreover, from the point of view of the energy, the 
functions Vj are truly lower order, i.e., they will satisfy 



I 



\V t ,xVj{t,x)\ 2 dx = 0{f ) t — (2.6) 
for all j > 1. In contrast, one of course has 

\V t/X u (t,x)\ 2 dx ~ 1 f->0 



J 



We shall now focus on the first two steps of the construction, i.e., u -- 
uq + v\ + V2- Let us compute the error resulting from Uq. Define D : = 
I + rd r = \ + Rd R . Then 



eo := ZquinticUo = ^ (t) [ ( j) % (f) (£> 2 W) (R) + (j)'(t)(»W)(R) 



f 2 e =:AJ(t) 



1 - R2/3 9 _ 3QR 2 + R 4 

tUl" 7— 7 + 



(2.7) 



(l + K 2 /3)5 (l + R 2 /3)5 

Here cuy are nonzero constants depending on v whose values do not concern 



us. 



2.3. The first correction. Then t 2 e = A{t)^0(R 2 (Ry 3 ) as R -> 00. This 
error blows up as f — > like f -2 . The goal is now to reduce it — in fact turn 
it into an error that vanishes as t — ► — by adding corrections to Uq, the 
first one being v\. We will do this by setting A 2 (t)LQV\ = e§ where 

L := d\ + |d R + 5W 4 (R) (2.8) 

Note that this is the linearized operator obtained by plugging «o + V\ 
into (|2.1|) and discarding d t altogether. While this may seem strange, the 
idea is to look first at the regime < r « t where dt should matter less 
than 8 r . We shall see shortly that V\ has the good property that it decays 
like (fA(f))~ 2 , but it produces errors for the nonlinear PDE that grow in r 
too strongly. To remove this growth, we carry out a correction at the sec- 
ond stage with a differential operator near the light cone r = t. There the 
self-similar variable a = ? becomes important. 

Now we discuss v\ in more detail. A fundamental system of Lq is 

, . l-R 2 /3 , s l-2R 2 + R 4 /9 

(pi (R) := '—, <p 2 (R) := '- (2.9) 

(l + R 2 /3)i R(l + R 2 /3)i 



The operator 



Lq = RL R- 1 = d\ + 5W 4 (R) (2.10) 



FULL RANGE OF BLOW UP EXPONENTS 5 

has a fundamental system 

(211) 

l-2R 2 +R 4 /9 ~ 

<p 2 (R) := R(p 2 (R) = f- = Rip 2 (R~ 2 ) 

(1 + R 2 /3)5 

The right-hand sides here are for large R, and the \pj are analytic around 0. 
The Wronskian is 

y' 1 (R)y 2 (R)-y 1 (R)f' 2 (R) = l (2.12) 
Define p(£) := tA(t), and 

^{^LoV! = t 2 e , Vl (0) = v[(0) = (2.13) 

We claim that 

P!(t,r) = n- z {t)L- x fe Q = A^(t)ii- 2 {t)0{R) as R -> oo (2.14) 
To be more specific, write 

t 2 e = A? (f) g! (R) + co 2 g 2 (R)) (2-15) 
see (|2.7|). Note that the g; are of the form 

gj (R) = R-'cPjiR- 2 ) R»l (2.16) 
where <pj is analytic around 0. ThenLo/y = with/y(0) = /'(0) = satisfies 

fj(R) = R- 1 J o ^(ROR'^XR') iR' - <p 2 (R) (f>i(R r )R r gj(R') dR' 

(2.17) 

for / = 1, 2. The one checks that 

fj(R) = hjR + b 2j + b 3j -^ + 0(1/R) as R - oo 
/y(R) = ci/jR 2 + 0(R 4 ) as R -» 

In fact, around R = the fj(R) are even analytic functions, whereas around 
R = oo one has the representation 

fj(R) = R(hj + b 2j R- 1 + R- 2 logR cp lj (R- 2 ) + R- 2 cp 2j (R- 1 )) 
=:R(F,(p) + p 2 G ; (p 2 ) log p) 

where <p\j,(p 2 j and Fj, Gj are analytic around zero, with p := R _1 . This 
follows from (|2.11[) , (12.161 ), and l|2.17|) . For future reference, we remark that 
the structure in (12.191 ) is preserved under application of D. In particular, 

v 1 (t,r) = A*(%- 2 (f)(<yi/i(K) + co 2 f 2 (R)) =: A*(f) p~ 2 (f) / (R) (2.20) 
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Define 

Ui := u + v t = (t) (W(R) + jT 2 (f) /(£)) 
In view of (|2.18[) , and R ^ /i (recall that we are inside of the light cone r ^ f) 

Ul (t,r) = A5(f)0(R- 1 ) R 1 (2.21) 

uniformly in < t < 1; moreover, we may apply fdf any number of times 
without affecting this asymptotic property. Finally, A(f) ~ 2 ui (£, r) is an even 
analytic function around R = 0. 

2.4. The error from «i. Setei := £.q U intk{ u i) ■ Then 
One has 

f 2 A~5(£)ei = A^5(f)((^) 2 - ^ f )(A5(f)wi(^rA(f))) - fi 2 (f)(10W 3 (R)o; 2 (f,R) 
+ 10W 2 (R)w 3 (£,R) + 5W(R)w*(t,R) + wfaR)) 

(2.23) 

where W\{t,R) = /i~ 2 (£)/(R). Then the nonlinearity in (12.231 1 is 

//(fXlOW^R^^R) + WW 2 (R)w 3 1 (t / R) + 5W(R)w*{t,R) + wfaR)) 
= Li- 2 (t)(lOW 3 (R)f(R) + 10W 2 (R)^i- 2 (t)f( y R) 
+ 5W(R)ii-\t)f\R) + Li- 6 (t)f(R)) 

(2.24) 

whereas 

A-5(f)((^ t ) 2 -^)(A5(0^i(£,rA(0) 

= ((fd t - (1 + v)£>) 2 - (td t - (1 + v)£>)) wi(£,R) 
Now [i(t) = £~ v whence 

f 2 (0 + ^1>) fT 2 (0/(K) = (2v - (1 + v)0)/(X) (2.26) 

Note that this is again of the form f(R) with / as in d2T8l) , d2T9l . Thus we 
can write 

?\-l(t)ei(t,r) = fT 2 (f)(/(R) - (10W 3 (R)/ 2 (R) + 10 f r 2 (f)W 2 (R)/ 3 (R) 

+ 5^\t)W(R)f(R) + ^(t)f(R))) 

(2.27) 
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We let a = j = jj = Rb, b := p -1 and isolate those terms in (|2.27|) which 
do not decay for large R. Since we are working inside of the light-cone, we 
have ^ a ^ 1. Now, abusing notation somewhat, 

ir 2 (t)f(R) = b 2 R(F(p) + p 2 G(p 2 ) log p) = ba{F{p) + p 2 G(p 2 ) log p) 

li- 2 (t)W 3 (R)f(R) = b 2 R- 3 n(p 2 )R 2 (F(p) + p 2 G(p 2 )logp) 2 

= b 2 R-\F(p) + p 2 F(p) logp + p 4 G(p 2 ) log 2 p) (2.28) 

p.-\t)W 2 {R)f{R) = b'R- 2 n(p 2 )R 3 (F(p) + p 2 G(p 2 )logp) 3 

= b 3 a(Hp) + P 2 HP) log p + p 4 F(p) log 2 p + p 6 G(p 2 ) log 3 p) 

where F, G can change from line to line. Similarly, 

p.- 6 (t)W(R)f(R) = b 6 R- 1 Q(p 2 )R 4 (F(p) + p 2 G(p 2 )logp) 4 

= b 3 a 3 (F(p) + p 2 F(p) logp + p 4 F(p) log 2 p 

+ p 6 F(p) log 3 p + p 8 G(p 2 ) log 4 p) 

p- 8 (0/ 5 (R) = fc 8 R 5 (F(p) + p 2 G(p 2 )logp) 5 

= b 3 a 5 (F(p) + p 2 F(p) logp + p 4 F(p) log 2 p 

+ p 6 F(p) log 3 p + p 8 F(p) log 4 p + p 10 G(p 2 ) log 5 p) 

(2.29) 

From (|2.28|) , l|2.29[) we extract the leading order 

t 2 \-^{t)e\{t,r) := /z -1 (t)(cia + c 2 b + (c 3 a + c 4 a 3 + c 5 a 5 )fc 2 ) (2.30) 

Indeed, from the first line in 02.281 ) we extract bfl(F(0) + pF'(O)) = bac\ + b 2 C2, 
whereas from the fifth we extract fc 3 aF(0). From the second line in (|2.29|) we 
retain fr 3 fl 3 F(0), and from the fifth one b 3 a 5 F(0). The point here is that with 



oo 



this choice of e x one obtains a decaying error as R 

£ 2 A-^)( ei -e;)(f,r) 

_ 7 rlogR 1 1 ( 13V) 

= V- W ( a ' b 'P lo S P'P) + x®2(a,b,plogp,p) 

where <&j(a,b,u,v) are polynomials in a, b and analytic in u, v near (0,0). 
Writing b = # we may delete the terms involving b 2 = ba/R on the right- 
hand side of fl2.30D , since they are of the form 02.311) . Thus, it suffices to 
consider the simpler leading error 

f 2 A-5(f) e 0(£, r ) := cxayr\t) + c 2 p~ 2 (£) = c x ab + c 2 b 2 (2.32) 

2.5. The second correction. Now we would like to solve the corrector prob- 
lem "near r = t", i.e., 

£{v tt -v n --v r )=-?e\ (2.33) 
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Note that we have discarded the nonlinearity on the left-hand side since it 
decays near r = t. This is designed exactly so as to remove the growth in R. 
We seek a solution in the form 

v(t,r) = A^Or^OfiCa) + ii- 2 {t)q2{a)) (2.34) 

with boundary conditions ^i(O) = 0,^(0) = and qi{G) = 0, q' 2 (0) = 0. 
These translate into the boundary conditions v(t, 0) = 0, 8 r v(t, 0) = 0. This 
will essentially be the function vi. In view of 

1 1 {a — h)v — i 

A(ty-^ a d t \(t)*ir a = d t + - — 2 - — 2 - 



we are reduced to the system 
a 



and 



Now 



t 2 (- {d t + ^^-f + Sn + - r d r )q z {a) = c 2 (2.36) 



t 2 (-{d t + ^) 2 + d rr + h)f(a) 
= ((1 - a 2 )d 2 + (2(j8 - 1)« + lxT x )d a - f + jS) f{a 



Define 



L^j := (1 - a 2 )d 2 a + (2(J3 - 1)« + 2a _1 )5 B - |3 2 + p (2.37) 
The system (12.351) , (12.361 1 therefore becomes 

Lv-iqi= da, L3v-iq 2 = c 2 (2.38) 

with boundary conditions ^i(0) = 0,^(0) = and ^2(0) = 0, q' 2 (0) = 0. 
Since for any integer n > 2 

Ljg a" = n(n + l)fl"~ 2 + (2n(jS - 1) - n(n - 1) + j3 - f)a n 

we see that (12.381) has power series solutions 

00 00 
Hi («) = 2 c i; fl2/+1 ' ^ 2 ( fl ) = 2 C2 /' ^ ( 2 - 39 ) 

;=1 ;=1 

convergent in \a\ < 1 and unique. In fact, C21 = ^ and c\\ = fy. 

Next, we determine the behavior of these functions as a — > 1—. We 
remark that in our case, always j6 > — \. There exists a fundamental system 
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of Lp of the form 

oo 

(2-40) 

i/> 2 («) = (1 - + J] & (1 - a/] = (1 - «)^+V°(a) 

f=i 

provided /S g Z^", see HI Lemma 3.6]. In that case these series define entire 
functions. On the other hand, if /3 e Z^, then i/^i is modified to 

oo 

jftoo = i + E ^ + c^ca) iog(i - «) 

with a unique choice of c. Notice that in this case the singularity of i/>i near 
a = 1 is no worse than (1 — a)P +1 log(l — a). 
In either case the Wronskian 

w = i/4 ip 2 - i// 2 

satisfies 

2(j3 - l)a + 2a- 1 . 

whence 

W(«) = fc(l - fl 2 )^<T 2 , fc * 
We define the Green function 

Gj}(a,a') := i/>i(a)i/> 2 («0 - \\>\{a')^ 2 {a) 
Then a particular solution of Lp q = f is given by 



q(a) = f Ga{a,a)W(fi)- 1 (l-&)- 1 f(a)dS 
Jo 

= jfc _1 f G«(fl,a)(l -a 2 y^ +1) a 2 f(a)da 
Jo 



fc-Vi(fl) ^2(2) (1 + a)-^- x #f(a)da 
Jo 

- r> 2 («) f - s 2 ) _tf+1) ?/(«) ^ 

Jo 



Returning to (12.381 ) we need to set /(a) = 1 and /(a) = a here. Note that 
q(a) as given by (|2.41[) satisfies the boundary conditions q(0) = q'(0) = 0, 
whence it agrees with q 2 and q\, respectively as given by (|2.39l) . 

Let us first assume that jS g Z^ . Then the term on the third line of (|2.41|l is 
analytic near a = 1. To analyze the behavior of the expression on the fourth 
line as a — » 1, we note that up to an analytic factor near a = 1 it equals 

(1 - af +1 f(l - a)-?- 1 h(a) da (2.42) 
Jo 
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where h is analytic in a neighborhood of [0, 1] . By inspection, (|2.42[) is of the 
form 

(1 - af +l ((1 - fl)^Hi(s) + c)H z {a) 

where c is a constant and Hi,H 2 are analytic near [0,1]. It follows that 
ft Zq, which means that v is neither an odd positive integer nor of the 
form ^±1 with n e Z+, one has that qi(a),q 2 (a) are analytic in the disk 
\a\ < 1 and near a = 1 they are of the form 

Qa(fl) + (1 - af +l Q 2 (a) (2.43) 

where Qj are analytic near « = 1. These functions have the property that 
after applying d a they remain in L 2 (0, 1) due to > —\. Evidently, they also 
become smoother as j6 increases, but they are never infinitely smooth (since 
|S is not an integer). 

On the other hand, if jl e Z^ , then the representation (|2.43|) needs to be 
modified with logarithmic factors log(l — a). We leave these details to the 
reader. 

Using a = Rfj.^ 1 we may rewrite (|2.34|) in the form 

vl{r,t):=^{Rq 1 {a) + q2 {a)) (2.44) 

where we have set q\{a) := a~ l qi{a). This ensures that both q\ and q 2 
have even power series in a around a = 0. Also note that these functions 
are 0(a 2 ) as a —*■ 0. We make one more adjustment: in (|2.44[) one has an 
odd expansion around R = 0, namely just the linear term R. We prefer to 
modify (|2.44[) as follows so as to retain the even expansion at R = 0: 



i 



v 2 (r,t) := ^-{-(R 2 (R)-%(a) + q 2 {a)) (2.45) 



A(*)i 

Note that for large R this captures the R growth of (|2.44|) , and the next 
order correction is R _1 . With this definition of v 2 we set u 2 := u\ + v 2 = 
Uq + V\ + v 2 . By construction, v 2 (t, r) is analytic in R around R = with an 
even expansion. Finally, (|2.21|) remains valid for u 2 as well. In other words, 
iio gives the main shape of the profile as a function of R. 

2.6. The error from u 2 . We define 

^2 := -Cquintk{u 2 ) = -Cquintic{Ul + V 2 ) 



2 

£quintic(Ul) + u\ - {U\ + V 2 f + (d tt ~ d„ - -d r )v 2 

e% — e® — 5u\v 2 — lQu\v\ — 10u\v^ — bu\v\ — v\ 



(2.46) 



2 



+ (dtt-drr- jd r )(v 2 -V° 2 ) 
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We determine f 2 A(f) ?e 2 . First, from (I2.31D 

t 2 A-5(0( ei - e ?)(^r) 

rlogK 1 i ( 2 - 47 ) 

= P W[^-*i(«/fe/plogp/P) + ^^(s, &,p log p,p) 

for R ^ 1. For |R| < 1 we read off from $227} and JZ321) that 

f 2 A-^)(e! - e ?)(t,r) = 0( M - 2 (t)) (2-48) 

This holds uniformly for small times, and tdt can be applied any number of 
times without changing this asymptotic behavior as R — > 0. 
Next, for large R 

t 2 A-^{t)u\v 2 = 0(R- 3 a 2 ) = OiR- 1 ^)- 2 ) 
which is of the form (|2.47|) . The final nonlinear term contributes 

t 2 A-^(t)v 5 2 = ^ 8 0(R 5 ) = ^R^Oi^R 6 ) = 0{R-^i(t)- 2 ) 

We leave it to the reader to verify that the other nonlinear terms behave in 
the same fashion. For small R, the nonlinear terms are 0(/i~ 2 (£)R 2 ). 
The difference V2 — contributes this: for R > 1, 

A(t)-?t 2 (d tt - d rr - h r )^((R - R 2 <R>- 1 )^ 1 («)) = 0(R-V(r 2 ) 

(2.49) 

Here we used that cj\(a) = 0(a 2 ) for small a. For small R this term 
is 0(p,~ 2 (£)R). By inspection, ei = -Cq U intk{ u 2) has an even analytic ex- 
pansion around R = 0, and by the preceding we gain a factor \i~ 2 for all R, 
and the decay is at least as R — > oo. 

2.7. Iterating the construction: the next corrections v 3 and V4. We now 

return to Section [231 in which we constructed V\ from eo- Here we need 
to determine ^3 from e 2 via the same route, i.e., by solving A(t) 2 LoV3 = e 2 
with zero Cauchy data at R = 0. The only essential difference between 
e$ and e 2 is that the latter looses a logR in terms of the R _1 decay. The 
dependence on a, b makes no difference in terms of the asymptotic behavior 
or the construction, however. Therefore, the same construction as before 
yields 

A(Q* 

as R — > go. Moreover, the asymptotic expansion involves terms with both 
log R and log 2 R. Around R = we have an even Taylor expansion starting 
off with R 2 . The 

u 3 := u Q + vi + v 2 + v 3 

still obeys the decay law (|2.21|) : in fact, any finite number of powers of log R 
can be lost since they are more than made up for by the gain of an extra 



v 3 (t,r) = -^-O(RlogR) 
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power of fi~ 2 in V3. In this sense, V3 is of strictly smaller order than both 
V\ and v 2 which can be comparable to uq on the light-cone r = t. For v\, 
one repeats the same construction that lead to v 2 above designed to obtain 
a decaying error. The error will be on the order of ,u~ 4 . 

The process can then be repeated any number of times. In [8 J the authors, 
together with D. Tataru, formalized this using various function algebras 
designed to hold the Vj the associated errors ey. The details are as follows. 
We have 

ttjfc = V k + Ufc_i 

The error at step k is 

2 

e k = (-df + d 2 + -d r )u k + u\ 

we define v k by 

(d 2 + ^d r + 5ulj v 2k+ i + 4 = ( 2 - 5 °) 

respectively 

+ d 2 r + h)j v 2k + e\ k _ x = (2.51) 

both equations having zero Cauchy data at r = 0. Here at each stage the 
error term e k is split into a principal part and a higher order term (to be 
made precise below), 

e k = e\ + el 

The successive errors are then computed as 

C2k = 4-1 + N 2icM> C2k+1 = 4fc _ 8 f V 2k+l + N 2k+ l(v 2k+1 ) 

where 

Na+i (v) = 5(u 4 2k - ul) v + Wu lk yl + Wu 2k v 3 + 5u 2k ^ + v 5 ( 2 -52) 
respectively 

N 2k (v) = 5u* kl v + Wu 3 2kl v 2 + Wul kl v 3 + 5m 2 h^ 4 + v 5 (2.53) 
The function spaces are as follows. First, the one relevant to the Lq iteration. 

Definition 2.1. S m (R k (logR)^) is the class of smooth functions v : [0, 00) — > R 
with the following properties: 

(i) v vanishes of order m and R~ m v has an even Taylor expansion at R = 0. 

(ii) v has a convergent expansion near R = 00 of the form 

00 

v(R) = J]J]c lj R k -'(\o S Ry 

The importance of even expansions in R near zero lies with the fact that 
only those correspond to smooth functions in R 3 . For the same reason, we 
will work with even m. Second, we introduce the space arising from the 
Sturm-Liouville problem near a = 1. 
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Definition 2.2. We define Q to be the algebra of continuous functions q : 
[0, 1] — > R with the following properties: 

(i) q is analytic in [0, 1) with an even expansion at and with ^(0) = 0. 

(ii) Near a = 1 we have an expansion of the form 

oo oo 

q{a) = q (a) + £(1 - af® +1 £ 9 y(«)(log(l - a))' 

with analytic coefficients qo, qif, if v is irrational, then q^ = if ; > 0. The 
/3(f) are of the form 

^ ((2* - 3/2)v - 1/2) + 2 ((2k - l/2)v - 1/2) (2 54) 

where X' are finite sets of positive integers. Moreover, only finitely many 
of the qij are nonzero. 

We remark that the exponents of 1 — a in the above series all exceed \ 
because of v > 0. For the errors we introduce 

Definition 2.3. Q' is the space of continuous functions q : [0, 1) — » K with the 
following properties: 

(i) q is analytic in [0, 1) with an even expansion at 0. 

(ii) Near a = 1 we have an expansion of the form 

co co 

q(a) = q (a) + £(1 - «)«'"> £ «fo-(a)(log(l - a))>' 
i=l ;'=0 

with analytic coefficients <^o, ^i;, of which only finitely many are nonzero. 
The jS(z') are as above. 

By construction, Q a Q'. The family Q' is obtained by applying ar l d a 
to the algebra <3. The exact number of log(l — a) factors can of course be 
determined, but is irrelevant for our purposes. 

Definition 2.4. a) S m (R k (log R) e ,Q) is the class of analytic functions v : 
[0,oo) x [0,1] x [0,b ] ^ E so that 

(i) v is analytic as a function of R, b, and 

v: [0,00) x [OA] ^Q. 

(ii) v vanishes of order m relative to R, and R~ m v has an even Taylor 
expansion at R = 0. 

(iii) v has a convergent expansion at R = co, 

co t+i 

v(R,a, b) = 2 2 C <V(«' b ) R^'(logR) ; ' 

i=0 ;'=0 

where the coefficients Cjj(-,b) e Q, and Cij(a,b) are analytic in b e [0,bo] for 
all < a < 1. 
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b) IS m (R k (log R) e ,Q) is the class of analytic functions w on the cone Co 
which can be represented as 

w(r,t) = v(R,a,b), v e S m (R k (log R) e ,Q) 

The same holds with <3' in place of Q. 

As in the first two steps we shall exploit that R, a, b are dependent vari- 
ables, by switching from one representation to another as needed. We shall 
prove by induction that the successive corrections and the corresponding 
error terms can be chosen with the following properties. There exist in- 
creasing sequences m^p^q^ of nonegative integers with m\ = p\ = 0,qi = 1 
so that for each k > 1, 

V2k-i e IS 2 (R (log R)'\ Q) (2.55) 

, i 

t 2 e 2k -i 6 IS (R (log R)K Of) (2.56) 

v * e is2 ( r3 ( lo s R ) ft 'Q) ( 2 - 57 ) 

e ^G^ 1 (logR)* Q) + b 2 IS°(R(logR)*Q')] (2.58) 

One can of course easily determine the optimal choice of m^, p^, cjk from the 
algorithm outlined below, but this is of no significance. We now inductively 
verify these claims. 

Step 0: The analysis atk = 
We observed in l|27)) that 

t 2 e e Ahs^R- 1 ) (2.59) 

as claimed. Now assume we know the above relations hold up to k — 1 with 
k 5= 1, and we show how to construct v^k-i, respectively v^, so that they 
hold for the index k. 

Step 1: Begin with e^^-r satisfying (|2.58D or 02.591) and choose i?2fc-l so 
J2551 froZds. 

If fc = 1, then define e° := eg- If > 1/ we use <2.58b to write 



where 



eik-2 = 4c-2 + e 2k-2 



mW 2 
KO 2 



(2.60) 



^^TTI^IS^RaogR)^,^) 
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We note that the term eL_ . can be included in e^-i/ cf- H2.56)) . We define 
»2ifc-i as i n d2.50D neglecting the a dependence of ejL._„. In other words, we 
choose to treat a as a parameter and the error resulting from this choice will 
then be incorporated in e^-i- 

Changing variables to R in (|2.50|) we need to solve the equation 

where the operator Lq is the one from above. Then (|2.55|l is a consequence 
of the following ODE lemma. We remark that the statement of Lemma I27T1 
is not optimal with respect to the number of logarithms in R. But for the 
sake of simplicity we choose this formulation. 

Lemma 2.1. The solution v to the equation 

L v = fe S^R- 1 (log R) e ), v(0) = i/(0) = 

with integer £ ^ has the regularity 

v e S 2 (R(logR)0 (2.61) 

Proof. This follows form the representation (|2.11|) , J2.16I ), and (|2.17|) . Indeed, 

v(R) = R- 1 (qk(R) £ q> 2 (R>)R'f(R') dR' - <p 2 (R) j ^(R')R'f(R') dR') 

J* oo e+i 

y?2(y- 2 )22c ;7 i/- ! '(logy)Mi/ 
Ro ;_n ;_ n 



(=0 ;=0 

oo e+i 

-<MR- 2 ) ^(y- 2 )2SQ;y^'(iogy)^y 

jR o !=0 ;=0 

+ ct^R) + c 2 M R ) 
The first integral in (|2.62|l contributes 

r-^r- 2 ) ^^^SS^y^'^sy) 7 ^ 

jR o n=0 ;'=0 /=0 

oo oo f+i 

= const ^(R) + i^i (R- 2 ) 222 Xli-m-it-i^nijR'-^ilogRy 

n=0 (=0 7=0 

oo oo t+i 

+ R-%(R- 2 ) 2 2 S *[i-2»-i=-« (logR)^' +1 



(2.62) 
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which lies in S 2 (R (logR)^), whereas the second one contributes 

f 8 oo oo t+i 

2 b «y~ 2n 22 c vy"( lo gy) / ^ 

jR o „=o !=0 j=o 

oo oo t+i 

= const cp 2 (R) + MR' 2 ) 222 X[-2n-i*-i]b„ ij R 1 - i - 2n (lo S R)i 

n=0 (=0 ;=0 

oo oo t+i 

+ u^ 2 ) 222 x[-zn-i=-i-\Kij (io g Ry +1 

which again lies in S 2 (R (log R-Y). □ 
We remark that V\ as constructed above satisfies 

see d2l8l) . 

Step 2: S/zoa; that ifvji-i is chosen as above then (|2.56|) holds. 
Thinking of i?2fc-l a s a function of t, R and a we can write e^t-l in the form 

ex-l = N 2k -i(v 2 k-i) + ESic-i + E a v 2k -i 

Here A^-i^ic-i) accounts for the contribution from the nonlinearity and 
is given by (|2.52D . E t v 2 h-\ contains the terms in 

- d tt V2k-i(t,R,a) (2.63) 

where no derivative applies to the variable a, while E"v2k_i contains those 
terms in 

(d ft ~ d„ - -d r ^V2k-i(t,R,a) 

where at least one derivative applies to the variable a (recall that in Step 1 
the parameter a was frozen). We begin with the terms in N^-i- We first 
note that, by summing the V t over 1 < j < 2k — 2, 

x\ 

u 2k - 2 -u e— IS 2 (R (log R) n , Q) (2.64) 



for some integer n = n(k) ^ 0. The first term in Nw-iivik-i) contributes 

f2 ( U 2)c-2 ~ U t) V 2k-l = t 2 [{u 2 k-2 ~ Wo) 4 + 4(«2ic-2 ~ "o) 3 "0 

-6(w 2 ic-2 - U ) 2 U 2 + 4(« 2 Jc-2 - W0)Wn] y 2lc-l (2-65) 
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Using (|2.64[) we compute 

£ 2 ("2fc-2 - «0) 4 f2)t-l 



l 



IS 8 (R 4 (logR) 4 "^)^ IS 2 (R (logR)'\<2) 



A 5 



« 6 IS 2 (R- 2 (logR) 4 "^)— IS 2 (R (logR)^Q) 



as well as 

a 



t (U 2 k-1 ~ Wo)"o y 2)t-l 

6 f 2 ^ IS 2 (R (log R)» Q)Ais°(R- 3 )^ IS 2 (R (logR)™* 
As 



The other two terms in 02.651 ) are similar. Next, compute 

^ e (^ ISl0(R5(1 ° SR)5m ' C ' Q) 

AlR6 IS 4 (R-i (logR) 5 '\£) 



(fA) ul 



(t\) 2k 

Aj 
(fA) 



fl 6 fo 8(fe-l) IS 2( R -1 (l og R) 5 "\<3) 

^(R- 1 (logR^G) 



and 



t 2 u\ k _ 2 v\ k _ x 



e A^(fA) 2 IS°(R- 3 ,Q)^ ¥ IS 4 (R 2 (logR) 2 "\G) 
13 6 a " 2 IS 4 (R- 1 (logR) 2M * / G) 



(fA) 



2fc 



A5 

' " is 2 (ir 1 (log a)**, Q) 



(fA) 



2Jfc 
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with similar statements for w2 jc-2 y 2Jc-i anc ^ Ulk ^ lV \k-Y Summing up we 
obtain 

This concludes the analysis of N^t-i (oafc-i)- We continue with the terms in 
E t V2k-i, where we can neglect the a dependence. Therefore, it suffices to 
compute 

f2d t ((^ IS2(R(1 ° gRr ' C) ) C (^ IS2(R(1 ° sK)mt) 
Finally, we consider the terms in E a v 2 k-i- With 

A5 



vzk-i(r,t) = -^w{R,a), w e S 2 (R (lo S R) m *,Q) 



we have 



-IRa^WRg - 2a~ l w a {R,a) + (a 2 - l)w aa (R,a) + 2aw a (R,a)] 
Since Q are even in a we conclude that 

(1 - a 2 )d aa , ad a , a- l d a : Q - Q' 

and therefore 

f^^e^IS^RaogR)^,^) 

This concludes the proof of ( 12.561 1. 

Step 3: Define v%k so that (|2.57|) holds. 
We begin the analysis with eafc-i replaced by its main asymptotic component 

at R = oo: 

= ^(E«?/(«)GogK)>'+& £ + «^(«)](logR)j 
P 7=0 y=o 

P)c+1 

+ fo 2 2 ?y(a)(logay), q u ff,%eGt 
7=0 

(2.66) 

This is chosen such that f 2 (e2ic-i — fik-i) consists of terms either decaying 
at least like R (logR) p * +1 , or else gaining at least a factor fo 2 . Consider the 
equation (|2.5H> with ^t_i on the right-hand side, 

t 2 (-d 2 + d 2 + h r )v 7k = -t Z f 7k -l 
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Homogeneity considerations suggest that we should look for a solution v& 
which has the form 

, i Pk Pk+t 

5* = H£i (£ <(«)(iogRV+b 2 2 w£(«)(io g Ry 

^ ;'=0 K=l,2 j=0 

j=0 

The one-dimensional equations for are obtained by matching the pow- 
ers of log R. This gives the system of equations 



+ * + 7*) 



= -2jt=rH"( fl )- F /( fl )) 

= ^(fl K - 1 ^(«)-Fj(«)),K = l / 2, 
r 



where 



Fy(fl) = (; + 1) [((v + l)v(2* - 1) + 2fl- 2 )W^ 1 + (a" 1 - (1 + v)a)d fl W^ 1 ] 



+ (/ + 2)(; + l)((v + l) 2 +«- 2 )wJ 2 

;'+i,k 



Fj(fl) = (; + 1) [(2fc(v + l)v + 20W£ + (a" 1 - (1 + v)fl)3 fl W; 
+ (; + 2)(/ + l)((v + l) 2 +fl- 2 )Wf 2 ' K 



>/+1,k 
2*: 



'2* 



Fy(a) = (; + !) 



((v + l)v(2* + 1) + 2a~ 2 )W J ^ 1 + (a" 1 - (1 + v)a)d a W, 



"ik 



+ (/ + 2)(; + l)((v + l) 2 +a- 2 )Wf 2 



'2* 



(2-67) 



Here we make the convention that WL = for j > p^ and WZ£ = WL = 
for j > pk + 1. Then we solve the equations in this system successively for 
decreasing values of ; from p^ to 0, respectively p^ + 1 to 0. 



20 JOACHIM KRIEGER, WILHELM SCHLAG 

Conjugating out the power of t we get 

e{ - [s, + {2k ~ 3/ f ~ V2 ) 2 + % + ^)w' t («) = «,;(») - W 

(2( _ ^ + (a - V 2)v- 1/2 )2 + ^ + ^ ^ _ ^ w _ JJw 



which we rewrite as an equation in the a variable, 

LpW[ k = aqjia) - Fy(a), fi = (2k- 3/2)v - 1/2 

L p W^ = a K - 1 ^(a)-F^a), = (2k - l/2)v - 1/2 ( 2 .68) 

LpW^ = %{a) - F>), f} = (2k + l/2)v - 1/2 
where the one-parameter family of operators is defined as above, i.e., 

Lp = (1 - a 2 H 2 + 2(fl- 1 + fljS - «)5 fl - ft + p (2.69) 

see d2.37| >. We claim that solving this system with Cauchy data at a = 
yields solutions which satisfy 



W'eaQ, j = 0,p k 



Wjea^Q, j = 0,p k + l,K = l,2, 
W ] 2k e <3, / = 0,p, f + l 



The choice of f> in H2.68D explains the appearance of the set 02.54I ). The fact 
that we need integer multiples of the exponents in (|2.68[) is a result of the 
nonlinearity which requires taking powers (moreover, Q is not an algebra 
otherwise). 

The claim H2.70[) is established as in the computation of vi above, see © 
Lemma 3.6] and [7, Lemma 3.9] for details. 

As in the case of vi, we need to make some adjustments. First, we modify 
V2k to ensure an even expansiorQ around R = 0: 

^=4( R2 < R r 1 2^ 1 <( fl )( 1 °s R ) i + E ^2>)( lo s R ) ; ' 



^(R)- 1 ^ a^wgWClogRy + 6 2 W^(«)(logR) 



y=o y=o 



We use the notation <R> = yl + R 2 . 
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This is not admissible because of the singularity of logR at R = 0. We 
therefore modify this expression further: 



V2k 

X 1 Pk 1 Vk+1 1 

= ^R\ R y^a- 1 W ] 2k (a)(-\o g (l+R 2 )y+ £ log(l + R 2 )) 1 

+ R 2 <R>- 1 2 a" x Wg (fl)(-log(l+R 2 )) ; + fc 2 W^( fl )(-log(l + R 2 )y) 

;=0 7=0 

Step 4: For v 2k defined as above we show that the corresponding error e 2k 
satisfies H2.58D . 



We can write e 2k in the form 



fe 2k = f (e a _i -e\_A + t 2 ie ^ 



-d 2 + d 2 + ^d r j v 2k J + t 2 N 2k (v 2 k) 



where N2/ c is defined by 02.531 ) and 

^4-1 - ^(2R 2 <R>^7(«)(^Og(l+K 2 )) ; 



7=0 
Pl+1 



2^(«)(-io g (i+R 2 )) ; 

7=0 

P/c+1 1 
7=0 

p,f+1 ~ 1 
^^(-logOL + R 2 ))'), 

7=0 



(2.71) 



with c/j,q/ ,qj(a) e C2'. We begin with the first term in e 2k , which has the 
form 

* 2 (^-i e ^[^(R-^logRF^ 2 ,^) +b 2 IS°(R(logRy\<2')] 

r 

The second term is admissible for (I2.58D . It remains to show that 

IS (R- 1 (log R)i*, a) c IS^R-^logR^O) + fc 2 IS°(R(logR)* / <3') (2.72) 
which can be seen as follows: for w e IS (R -1 (log R)*, <3') we write 

w = (1 - fl 2 )w + — - -R 2 w 
(fA) 2 
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Then 

(1 - a 2 ) w e IS (R- 1 (log R)#, Q), t^R 2 w e b 2 IS (R(log R) qk , €t) 

as desired. The second term in e 2 fc would equal if we were to replace 
\ log(l + R 2 ) by logR in both and z^; in addition, as in the case of ei 
above, error terms arise due to the replacement of R with R 2 (R) _1 . We 
leave those latter terms to the reader. To illustrate the former, consider the 
difference which is obtained upon replacing the derivatives of \ log(l + R 2 ) 
by derivatives of log R in the expression 

t2 + d ' + \ d ) f 2 W 2,(«) Ml + **) 

Computing these differences one finds that the second term in e 2 / c is a sum 
of expressions of the form 

— TSE-T 2 [s(R- 2 )(log(l + ^ 2 )) / - 1 + S(R- 2 )(log(l + R 2 ))'- 2 



(fA) 



S(R- 2 )(log(l + R 2 )y- 1 



a 

Since are cubic at it follows that we can pull out an a factor and see 
that this part of the error is in 

-i-IS^R-HlogR)^,^) 

which is admissible by (|2.72[) . The nonlinear expression is left to the reader. 

In summary we arrive at the following result. 

Proposition 2.2. For any N ^ 1 there exists a positive integer k such that u k as 
constructed above satisfies 

f \(£ q uinticU k )(r,t)\ 2 r 2 dr = 0(t N ) f -> 0+ 

1 

Moreover, u k - W A( - f ) e ^-O(R), wzf/z R = rA. 

The following proposition is going to be the key for the proof of Theo- 
rem ll.il It states that the approximate solutions of the previous proposition 
can be turned into actual solutions inside the light cone. 

Proposition 2.3. Fix N > 1 sufficiently large. Then there is k ^ 1 sufficiently 
large and there exists a function e(t, r) with 
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such that the function 

u(t,r) = U2k-i(t,r) + e(t,r) 
is a solution of inside the light cone. 

Proof of Theorem \T?l\ assuming Proposition 12.31 As in [8] one observes that u k 
can be extended outside the light cone such that 



f [\V x u k \ 2 + |V,W A(f) | 2 + (d t W A{t) ) 2 + \V t u k \ 2 ]dx < 5 

Jr>t 



for any prescribed 5 > 0, provided < f < fo with to = to(5, k, v) sufficiently 
small. But then, picking to sufficiently small, we can arrange that the 
H 1 x L 2 -norm of («(fo, -),u t (to, •)) is less than 6 in the region r > to- Then 
Huyghen's principle and the small energy global regularity imply that the 
corresponding solution remains of class H 1+ 2~ on the exterior of the light 
cone. □ 

The remainder of the paper is devoted to proving Proposition 12.31 The 
idea is that up to smoother errors, the principal source term ^-l ( we take 
u 2k-i t° be the approximate solution) may be reduced to an explicit algebraic 
expression. In fact, due to Huyghen's principle, we can and shall modify 
e 2k-i outside the light cone to simplify the analysis a bit. We shall then see 
that the error term e(t, r), which will be constructed via a suitable iteration 
scheme, can be written as a sum of a smoother and a "small" term (in the 
sense of amplitude). This decomposition in fact results naturally from the 
structure of the first iterate of the scheme used to construct e(t, r). 

3. Setting up the iteration scheme; formulation on the Fourier side 
Introduce the variables t = v~ x t~ v , R = A(f)r, and write 

e(T,R) :=Rf(f(T),r(T,R)) 
Then we get the following^ equation for e: 

(<9 T + AA^RdR^e - {d T + AA^Rd^e + £e 
= A- 2 (T)R[N 2k ^(e) + ea.!] + d T (AA^ 1 )e; ^ 
where the operator £, is given by 

£ = -3l-5W 4 (R) 

and we have 

RNa-iO) = 5(i4_! - u*)e + RN(u2k-i, e), 
RN(u2k-i,e) = R(u 2k -i + -|) 5 - Ru 5 2k l - 5u\ k _ x e 



2 In fact, the term d T (AA ')f was inadvertently omitted in (8), which has no bearing on 
the proof there, however. 
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Introducing the operator, with j6 v (t) = j^y, 

D = a T + i 3 v (T)(R5 R -l), 
we can also write the above equation as 

&£+f} v (z)De + £? = A- 2 (T)[5(u 4 2k _ 1 -u*)e + RN(u 2k _ 1 ,^+Re 2k _ 1 ] (3.2) 

In the following, we shall freely borrow the facts about the spectral theory 
and associated distorted Fourier transform contained in J8] as well as |2|. In 
particular, we recall that there exists a Fourier basis (p(R, £) and associated 
spectral measure p(£) satisfying the asymptotic expansions and growth 
conditions explained in [8, Section 4] such that 

rCO 

?(t, R) = x d (T)(f) d (R) + x(t, Z)<t>(R, $)p(S) dt; 

Jo 

For the asymptotic behavior of <p(R, £) in various regimes see (|4.2|l and (|4.4|l . 
Here the functions x(t, £,) are the (distorted) Fourier coefficients associated 
with e, and (pd{R) is the unique ground state with associated negative 
eigenvalue for the operator £. We also note the important asymptotic 
estimates 

p(£) rK 4 « 1, p(£) £*, £ » 1. (3.3) 

as well as the fact that near £ = as well as £ = oo the spectral measure 
behaves like a symbol upon differentiation. We shall henceforth write 

Then proceeding as in [2J, in particular section 3.5 in loc. cit. which uses a 
variation on the procedure in [8|, we derive the following transport equation 

for x(t,£): 

(£>* +/3 v (t)£) t + |)x(t,£) = K(t,x) + /(t,£), (3.4) 
where we have 

= (-4 i Sv(T)ra T x- j g2( T )(7C 2 + [^,7C]+7C+^ 2 7C)x)^) (3.5) 

withjS F (x) = ^y, and 

/(t,£) = r(A- 2 (T)[5(4_ a -«*)£ + RN(« 2fc _ 1/ £) + Re 2fc _ 1 ])(£) (3.6) 
Also the key operator 

D T = d T +p v (x)2i, ^ 

and we have 

91 2£rl ( 5 4- ^ 
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Finally, we observe that the "transference operator" % is given by the 
following type of expression 

< 3 - 7 > 

with mapping properties specified later on. 

In the following, we shall take advantage of the observation that the 
equation 

{p\ + j3 v (t)D t + !)x(t, Q = /(T, Q = ( ^ W } ) 

can be solved completely explicitly; in particular, imposing vanishing bound- 
ary data at t = oo, we obtain the following expression for the continuous 
part x(t, £): 

Jt A3 (a) pi(£) L Jt j A 2 (a) 

(3.8) 

We shall justify this below. On the other hand, one immediately obtains the 
elementary implicit relation 

"CO 1.1 

z-a\ 

(3.9) 

f d (a) = f d (o) - p v (o)d T x d (o) 
To derive (13.81) , define an operator 



f ~ 1 1 1 

x,(t) = J H d (x,a)f d (o)da, H d (x,a) = ——\S, d \~^ e 2 



(Mf)(x,t) :=A-i(T)p(«)»/(^€) 
and £)i := 5 T - 2j3 v (T)£<9 { . Then 

M _1 £>iM = £> := a T + |Sv(t)51 c 

Second, define 

(S/)(t,5):=/(t,A- 2 (t)£) 
Then S _1 5 T S = D\ whence 

A^T 1 S _1 5 T SM = D 

It follows that 

(SM)-\d 2 T + p v (x)d T + A- 2 (x)£)SM = D Z + p v {x)D + £ (3.10) 
Finally, one checks that H(x, E) : = s\rv{E,^oj{x)), co(x) := $ T A -1 satisfies 

(d 2 T + fi v (x)d T + a- 2 (t)<^)h(t, z) = o 

Hence, 

H(£, a, x) = £~ 3 sin J A _1 (u)dtt), cj t 
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is the fundamental solution of the operator in the parentheses on the left- 
hand side of (|3.10[) . In order to measure the size of x, we use the norms 

= \\(0 a x\\ L 2 
dp dp 

The following lemma shall be used throughout the sequel without further 
mention: 

Lemma 3.1. Assume N is sufficiently large and write 

1/11^:= sup T^d/^OI^ + l/iCT)!) 
— dp T>TQ dp 

Then defining x(x, E) via (|3.8|) and x d {%) as the unique fixed point of (|3.9|) , we have 

\\ X \\ 2,*+\;N-2 ^ WfWj}^ 
L dp ~ dp 

In fact, for the discrete spectral part we can improve this to 

Mt)|<t- n ||/|| l 2,„ 

— dp 

The proof is a straightforward consequence of the fact that 



£ 2 sin 



< T 



2 ' 2 ~ 2 J ~ 

such that (I3.4I) admits a fixed point x(t, ■ ) e t~ n L ' 2 ,D t x(t, •) e t _n_1 L ' 2 



r i f ff i 
A(t)<^J \~\u)du 

What we shall effectively prove now is the following 

Proposition 3.2. Given N sufficiently large, < v < \, there is k sufficiently 
large and there exists a function ejk-i £ such that e 2 k-\\ r <t = ezk-l\r<t and 

' 2 

Proof of Proposition 12.31 assuming Proposition 13.21 This is now a consequence 
of the fact that T is an isometry from Hj^j to L^ a . The discrepancy between 
2a on the one hand, and a on the other hand is due to the frequency being 
£5 rather than £). □ 

Proof of Proposition 13.21 We shall solve (13.41 ) via an iterative scheme, namely 

{p\ + ]6 v (t)£> t + .(t, 5) = K(t,x m ) + /;-i(t,|), (3.11) 

with 

■&(-\(<l/L . , ^ , Pv. 



= ~4|3 v (t)'7C2) t x / _ 1 -^(t)^ + [^, , 7C]+7C+^'7CJx._ 1 (3.12) 

Pv 

and 

/;-i(t,£) = r(A- 2 (T)[5(4_j - «*)£)_! + RN(u2k-i,£j-i) + Rhk-i]) {£) 

(3.13) 

and we of course set 

ej(T,R) = x dtj (z)(f) d (R) + Xj(T,Z)<t>(R,Z)p(Z)dZ, j > 1, 

Jo 
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while we also set Eq = x Q = 0. In particular, /o = e F(Re 2 k-i), ^(T/*o) = 0- 
This underlines the importance of the first iterate e\ since it will determine 
the smoothness of subsequent iterates due to the smoothing properties of 
the parametrix. We next turn to a careful analysis of the first iterate. 

4. The first iterate 

i 

In light of the fact that f 2 e 2 jc-i 6 j§fk IS ° ( Rl+ > &) and the 

precise definition 

of this space from before, it is clear that e2fc-l is C°°-smooth except at the 
light cone r = t. Furthermore, by subtracting functions of regularity at least 
Hj R and choosing k large enough, we may replace e 2 fc— l by an expression of 
the forrr@ 

t 2 e 2k -i = C(T)T- N - 2 (log(l - a)y(l - a)^, C(x) = 0(1), (4.1) 
for a finite collection of indices i. More precisely, we can write 

t 2 e 2k -i = 2 Q(T)T- N - 2 (log(l - a))\l - «)¥ + Eat-! 

i 

with Ea-i e t- n H^ r . 

The principal issue now becomes what the effect of the parametrix (|3.8D 
on this kind of expression is. Then the key to proving Proposition |3.2| is the 
following result. 

Lemma 4.1. Define e^k-i to be the function obtained from 

CT- N - 2 (io g (i-fl)y(i-«) v -^u 

by truncation tor ^ f. Then the function x{%, £,) defined via 03.81 ) with 

f(T,Z) = \(z)- 2 T(Re 2 k-i(T,m) 
satisfies x = x^ + x^ 2 \ where 

•) e t- n L^-, ,( 2 ) e r- N L^, ^^(x®) e 
provided k is sufficiently large. Finally, we also get the bound 

D.xex-^L 2 ^, \x d (z)\ZT- N -' 
where x^ is defined via 03.91) , with f replaced by Re 2 k-i- 
Corollary 4.2. Modify e^-i by restricting the terms 

Q(T)T- N - 2 (log(l-fl))''(l-fl)^ 

to r < t while smoothly truncating Bjk-x to a dilate of the light cone. Then the 
same bounds as in the preceding lemma apply to x, the true first iterate. 



'Recall that we may modify e^-i arbitrarily outside the light cone. 
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Proof of Corollary \4~2\ This follows from Lemma l4~T1 as well as Lemma I3T1 
applied to E^t-i- □ 

Proof of Lemma I47TI The proof of the bounds for small £ proceeds exactly as 
in [8], and so we shall now focus on E, » 1. The expression fi^e^k-x ( T / ')) (£) 
is given by 

poo 

T(Re 2 k-x(?, ■))(£) = (j)(R,Z)Rezk-x(*,R)dR 
Jo 

and here we can restrict the integration to < R < vx. We shall next use the 
precise asymptotic expansion for (p(R, £,). Specifically for R£2 > 1, we get 



<K2U) = + «(£)/-(*,£), /-(R,4) = A(R,£) (4-2) 

where we have the expansions in Hankel functions 

f+(R,Z) o(R&,Z), 

1=0 

and ipt uniformly bounded in R with ip~j~(R) = 0(Rl). In particular, we get 

/ + (R / 0=« fl ^ i +O(H) 

for large E, and R£5 ^ 1- 

Another key ingredient is the precise asymptotic formula for the coefficients 
«(£). In fact, from 0, Cor. 3.7, we have 

«(0 = ^ r (i + o(H)), £»i 

Neglecting for now all terms of the form 0(E~^), we are then lead to the 
following integral: 

T(Re 2 k-x{c,-m) 

r-l rva i r-l rva i 

= ^-J e lR ^Re 2k -x(o / R)dR-^-j e^ 7 Re 2 k-x(o, R) dR 
i . 

= ^f^J o e^ R ^ va ^ J Re 2 k-i (a, R) dR 



2i 



f a e-' (R - va) ^R? 2 ,c_i(a,R)dR 
Jo 



Recalling the definition of e 2 / c _i and replacing the outer factor R by vo 
up to an error of type Eat-i, the preceding is seen to be equal to a linear 
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combination of terms of the form 

(logx)* sin(£5x 



m m 11 C va log xV sin 45* 

(loga) N2 a- Nl cos(va4 5 )^" 5 \ dx 

Jo IT 

= (loga) N2 a- Nl 4 _3 Tcos(va45)C(4) 

Nt>N + l, 

mm ii f va (logxV cos(£5x) 

(loga) N2 a- Nl sin(va45)4" 5 ; <fa 

Jo i- 

= (loga) N2 a" Nl 4 _ ^sin(va45)C(4) 
Ni > N + 1, 

for some uniformly bounded functions C(4),Cy(4) with symbol behavior. 
Note that the errors 0(o~ Nl ~~ <; _5 + ) are in a~^~ 2 L^ , which corresponds 

to H r2 as for the error terms E^k-i introduced further above. Now apply 
the parametrix H3.8D to the first two terms in the above two expansions, 
multiplied by A(a)~ 2 . Observe that 

A(t)£5 f \- l (u)du=VT(l-(-) 1 v)ES 

Jt O 

and so we obtain the respective contributions 

Jt A°+(a) A 2 (a) p i(,r) a 

T i+v- 1 

(log a) N2 a" Nl cos K (v —Z~ 2 ) da, k = 1,2, 



as well as 



Jt A°+(a) A 2 (a) p5 (^) a 

A+v- 1 



(loga) N2 a" Nl sin K (v— —^) da, k = 1,2, 
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where now Ni > N + 1 + 2(*± 1 ). Expanding the sine function in the 
integrand, we encounter in the case k = 1 the problematic terms 

~°°A 0+ (t) A 2 (t) P'(l^^)(loga) N 2 ^i+v- 1 

5 

cos(vt£< 



Sm(VT52)5 4+ J T A^) C( A^)^^ ^- C ° S ^ 



where the bacO terms cos 2 ^ 1 ^ ^^sin 2 ^ 1 ^ £2) may be replaced up 
to an oscillating term by |, in which case the preceding further simplifies to 

V)( n w thr-A: r^WOoga)^ A 2 (t) P^(||f) 

J T A u+ (a) a N i A 2 (a) p i(g) 

(4.3) 

as well as a similar term with cos(vt^z) in front. Here we have not gained 
the necessary decay for the Fourier coefficient x. On the other hand, this 

term has a well-defined oscillatory behavior in terms of £2. Thus for R < y, 
we get 

rCO rCO 

( p{R,Z)xM{T,t)p{Z)dt=\ X a <K1U)X (2) (T,£)P(£)# 
JO Jo L-Ke 2 S1J 

poo 



In the first integral we have the implicit phase functions e ±!VT ^ 2 , and per- 
forming integration by parts with respect to £2 leads to a gain of t -1 ^?, 

which makes the integrand absolutely integrable with respect to d^. In the 

1 

second integral one has the implicit phases e'(± R ± VT )5 2 / and here we gain 
(+R + vt) _1 ^ _ 3. More precisely, for the first integral one writes (see HJ) 

00 

= ^o(R) + R- 1 2(K 2 ^)^(R 2 ), (4.4) 
7=1 

with (pj{R 2 ) < R 2 , and if this expression is hit by a 5 1 , it changes to 

00 

2R^£;(R 2 £)/-> ; (R 2 ) / 



4 in the sense of non-oscillatory 
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and we have R 2 < RE,~? on the support of the integrand. The extra factor R 
is used to absorb the R _1 in 



The case when the operator d 1 hits the other terms in the integral are 
handled similarly. 

We have now reduced matters to controlling integrals of the forrrQ 

sin(vT^) f° A°+(t) ^A 2 (t) c . P^(g|fjO - T^\i, (logoff 
cos(vt,S5) f co A°+(T) A 2 (t) sHjMZ) t i+v-! (loga) N2 j 



Jt A°+(a) v A2(a) w p | ( £) 

with S = sin, cos and K e {1,2,3}. Here we perform integration by parts 
with respect to a, which gives the desired gain placing the term into 

dp 

To conclude the bounds for the first iterate, we still need to bound the 
contribution from all the errors which arose when we replaced (p{R,E) by 

cIm(z'<!;~2e !R £ 2 ), as well as the errors of type E^-i- Note that the former 
type of error contributes to the Fourier transform of Re-&_\ a term of the 
schematic form 

rva 

r'Aia)- 2 Re^ dR 
Jo 

and we can use a crude L^-bound for the integrand to infer that this 
contribution is bounded by < T~ N £,~i, which places this contribution into 

2 - — 

z~ N L d '* . The contribution of errors of type E^-i was handled in the proof 
of Corollary 14.21 This completes the proof of the lemma. □ 

5. The second iterate 

We observe that Lemma l4Tl implies the following: the first iterate E\ can 
be written as e% = + where 



1] hV ' 11 R Iil 2r 



I^DI] 3 + ||f_|| „< T -N 



dR 



(5.1) 

r ! -'\\ . . 4- !lv.„ .,i^L w < r~ s 

dR 



^These arise upon replacing sin 2 ^ 1 '*^ g?) by sin 2 (v Tl ^l 1 e\) — \ and similarly for 





ll ?(2) II 


T CO H" 


■ II R II 
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where M > 2 can be chosen arbitrarily large (with implicit constant de- 
pending on M). By radiality we then obtain 

To control the second iterate write 

+ i S v (T)D T +|)(x 2 -x 1 ))(t / ^) = ^(t,^) + Ai/xCx^) 

with 

Ax/iCa) := !F(A _2 (T)[5(«2 Jt _ 1 - «g)ei + RN(u 2fc _ 1/ £i)]) (£) 
Then we claim 

Lemma 5.1. We have the estimates 

(x 2 - *i)(t, •) e T- N L 2 / p , \(x 2 - Xl ) d (T)\ < z- N -\ 

D T (x 2 -x 1 )(z r ) e t-"- 1 ^, \d x (x2-xi) d (x)\ < t""" 1 

hold. In fact, one gains a factor ^ in the corresponding norm bounds. 

Proof of Lemma |5J1 From Lemma l3?tl we conclude that, with x, f as in (|3.8|) , 
we have 

SUpT N ||x(T / -)|l 2 , a+ l < SUpT N+2 ||/(T,-)|| L 2, 1I 

T>0 L T>0 

SUpT N ||D T x(T, •)|| L 2, a ^ SUpT N+1 ||/(T / -)|| L 2,« 

T>0 T>0 

Further, from (8), 13 we have the operator bounds (here a is arbitrary) 

K dc :L 2 d '^R, <K dc :R~>L 2 d >«, K dd : R - R 

In fact, the operator "7C C(; is a smoothing operator, and corresponds to the 
operator TCq in HI • Our strategy is to exploit the smoothing effect of the 
wave parametrix d3.8b - This indeed leads to a derivative gain for all the 
terms in ^?(t,x 1 ): assuming the that the Fourier transform X-, of the first 
iterate can be decomposed into two terms with bounds as in Lemma 14.11 
we get 

K{ % x x ) e T- N - 2 L 2 d j 
Application of (13.81) leads to expressions in t^ n L 2 '^ ; moreover, applying 

2 - 

D T to these terms leads to expressions in t _n_1 L^ 2 , as required. The 

contribution to the discrete spectral part is also immediate. 

We next turn to the contribution of Ai/i (t, E) . In fact, we claim that all these 
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terms lead to a contribution in L^' 1 upon application of (|3.8|) . To see this, it 
suffices to check that 

T{A- 2 {t)[5{u^_ x -4)8, +2W(« afc _i,£i)])(*) e t- n - 2 L 2 /> 

For the first term, use the relation from [8] that 

u^-i - uo = 0(— ) 
from which we easily infer 

|-T(A- 2 (t)[5(4_ i _ u 4 )?l]) |^ i < (Atr 2 ||?i|| H j R - ^ 2 Fill^ R 

dp 

which is bounded from Lemma l4~Tl To control the source term RN(«2jc-i, ?i ), 
we consider the two extreme possibilities 

70. 

U 2k-1R> (r) £ 1- 

To bound the first term on the left, we use 



J dR 



while for the second term we have 



' £ lyt~ II < n £ i ii5 , N £ i 1. 4 .I ~ || 

i< ""rfK i< L dR K L dR n dR 



and we have 



N-2 



Note that the rapid decay rate of £1 gives much more than t -decay. 

6. The higher iterates 

Here we repeat the procedure of the preceding section, except that we 
replace x_ 2 ~-\ ^Y*k~ -k-i'^ ^ 3,and we replace one copy of £iby £fc_i — £k-2 
in each source term. Then we can literally repeat what we did before. To be 
specific, we claim the following: 

Lemma 6.1. We have the bounds 

(x k - x^Xt, •) e T- N Lf p , \(x k - x k _ t ) d (x)\ < z- N -\ 

G T (x k - ^-i)(t, •) e t-^L 2 /; , \d T (x k - x fc _O d (T)| < t- n -' 
In fact, one gains a factor in the corresponding norm bounds. 
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The proof is by induction on k, and is in all respects identical to the one of 
the preceding lemma, except that there is no one factor e^-\ — £k-2 involved 
in the analogue of Ai/i(t, £). The factor (jj) k comes from the repeated time 
integrations. □ 

The fixed point of (|3.4|) is now found by iteration, which completes the 
proof of Proposition |2.3l □ 
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